Relativistic and Radiative Corrections to the Mollow Spectrum 
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The incoherent, inelastic part of the resonance fluorescence spectrum of a laser-driven atom is known as 
the Mollow spectrum [B. R. Mollow, Phys. Rev. 188, 1969 (1969)]. Starting from this level of description, 
we discuss theoretical foundations of high-precision spectroscopy using the resonance fluorescence light of 
strongly laser-driven atoms. Specifically, we evaluate the leading relativistic and radiative corrections to the 
Mollow spectrum, up to the relative orders of {ZaY and a [Zaf' , respectively, and Bloch-Siegert shifts as 
well as stimulated radiative corrections involving off-resonant virtual states. Complete results are provided for 
the hydrogen \S-2Pi/2 and IS-lPz/^ transitions; these include all relevant correction terms up to the specified 
order of approximation and could directly be compared to experimental data. As an application, the outcome 
of such experiments would allow for a sensitive test of the validity of the dressed-state basis as the natural 
description of the combined atom-laser system. 

PACS numbers: 12.20.Ds, 31.30.Jv, 06.20.Jr, 31.15.-p 



I. INTRODUCTION 

Experimental possibilities in high-precision spectroscopy 
have received a rather significant 'boost' in recent years due to 
the availability of phase coherent regularly spaced frequency 
combs that may bridge large frequency intervals between fre- 
quency standards and optical transition frequencies |1]. In 
general terms, highly accurate spectroscopy may lead to an 
experimental verification of known theoretical models of the 
physical process under study. Precision measurements — in 
combination with theory — allow to obtain accurate values for 
physical parameters or fundamental constants [2J. With in- 
creasing accuracy, one may even ask whether the so-called 
constants are in fact constant HHHSIHllHll. These 
results may as well be used as input to more applied physics 
as the creation of unit standards e.g. for time and mass. 

On the theoretical side, quantum electrodynamics is one of 
the most accurate theories known so far. In many previous 
studies, the S'-matrix formalism has been used to obtain accu- 
rate predictions for experimental investigations. The 5-matrix 
relates the distant past to the distant future without referring 
to the dynamics of the intermediate times and leads effectively 
to a somewhat static description 1 1 1 ] . 

Thus, the S'-matrix formalism cannot, a priori, lead to a sat- 
isfactory description of quantum electrodynamic corrections 
to dynamical processes, and the subject of this paper is to 
provide a first step in the direction of a high-precision theo- 
retical description of dynamical processes including radiative 
corrections, using a laser-driven atom as a paradigmatic ex- 
ample. Obviously, the treatment of radiative corrections to a 
dynamically driven atomic transition requires input from two 
different areas, which are laser physics and quantum electro- 
dynamics. While the two areas are related, there are a couple 



of subtle points to consider when a unified understanding of 
a specific problem is sought, whose nature inevitably requires 
concepts introduced within the context of either of the two 
areas. In particular, it is known that the description of dy- 
namical processes requires considerable care in the treatment 
of the gauge-dependence of amplitudes, and with regard to the 
physical interpretation of the wave functions used in the math- 
ematical description 

A classic textbook example for a dynamical atom-laser sys- 
tem, well-known in theoretical quantum optics j^^l, consists 
of the Jaynes-Cummings model of an atom that contains two 
relevant energy levels interacting with a single monochro- 
matic laser field mode |23]. Due to the driving of the laser 
field, the atomic population undergoes Rabi oscillations. The 
population is driven periodically from the upper to the lower 
state and vice versa. The emission spectrum of this process 
with a strong driving field is known as the Mollow spec- 
trum i24ll . This case of strong driving may easily be inter- 
preted in terms of the so-called dressed states. Laser-dressed 
states are defined as the eigenstates of the combined system 
of atom and driving laser field | 25] and have proven to be use- 
ful in countless cases of both theory and experiment, one of 
which is the Autler-Townes splitting |26]. 

When evaluating radiative corrections to the Mollow spec- 
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FIG. 1 : Diagrammatic representation of a radiative self-energy cor- 
rection to the laser-dressed atomic state. The double line corresponds 
to the electron bound by the nuclear Coulomb field. The jagged line 
denotes the additional dressing of the bound electron by the (strong) 
laser field. The self-energy of a laser-dressed. Coulomb-bound elec- 
tron is a quantum-field theoretic problem in the presence of two clas- 
sical background fields. 
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trum, it is natural to start from the dressed-state basis, which 
consists of the natural eigenstates of the (strongly) coupled 
atom-laser system rather than the bare atomic states. It might 
be assumed that in order to fully treat the Lamb shift of laser- 
dressed states, it would be sufficient to simply correct the 
energies of the bare states that enter into the formulas for 
the generalized Rabi frequencies by the "bare-state" (i.e., the 
usual, ordinary) Lamb shift. Indeed, the first investigations 
on the problem |27| revealed corrections to the dressed-state 
"quasi-energies" consistent with this assumption. However, 
recently, it was found that at nonvanishing detuning and Rabi 
frequency, the Lamb shift of dressed states is nontrivially dif- 
ferent from the bare Lamb shift 1 28. .29] . Thus the distinction 
between evaluations in terms of the bare- and the dressed state 
basis in fact has to be made. In the limit of vanishing detun- 
ing, the coincidence of the bare- and dressed Lamb shift ef- 
fects on the detuning is obtained only after a summation of a 
specific series whose leading correction term may be obtained 
by carrying out the calculation to second order in the atom- 
field interaction. A diagrammatic representation of a radiative 
correction to the dressed state is shown in Fig.^ 

Thus in this article we present a detailed and complete the- 
oretical analysis of the leading nonrelativistic and relativistic 
corrections to the Mollow spectrum, up to the relative orders 
of (Za)'^ and a (Za)'^, respectively, and of Bloch-Siegert 
shifts as well as stimulated radiative corrections involving off- 
resonant virtual states, and laser-field configuration depen- 
dent corrections. The purpose is to enable a direct compari- 
son between theory and experiment for a high-precision spec- 
troscopic investigation involving laser-dressed states. Such a 
comparison to experimental data would allow to address ques- 
tions related to the physical reality of the dressed states (and 
their "quasi-energy") on the one hand, and of the nature and 
the interpretation of the various radiative corrections on the 
other hand. As a promising candidate for the experiment, 
we identify the resonance fluorescence spectrum of a strongly 
driven hydrogen 1 S-2P transition, which to lowest order may 
be described by the standard Mollow spectrum. A coherent 
Lyman-a source i30l|3TI has recently become available as a 
driving field, and we show that ionization into the continuum 
does not prohibit an experimental implementation. In partic- 
ular, we discuss corrections which are due to resonant and 
off -resonant excitations as well as the Bloch-Siegert shift, and 
corrections to the transition dipole moment and to the secular 
approximation leading to the Mollow spectrum. As a result, 
we provide theoretical predictions which are directly compa- 
rable to possible experimental data. 

The article is organized as follows. In Sec. HI] we intro- 
duce our system of interest and provide the relevant theoretical 
background for the further analysis. In Sec.|in| we evaluate 
the corrections to the Mollow spectrum, which we divide into 
modifications of the detuning (Sec. IIII A> and the Rabi fre- 
quency (Sec. nil St . The dominant relativistic corrections are 
of the order of {Za)"^, where Z is the nuclear charge number, 
and leading radiative effects lead to correction terms of the 
order of a (Za)'^ \n[{Za)^'^]. In Sec. lIVI we provide numer- 
ical data for the hydrogen lS-2Pj transitions (j = 1/2,3/2). 
Sec.|V]discusses and summarizes the results. 



II. MOLLOW SPECTRUM 

In this section, we introduce our system and recall results 
of previous studies which will serve as the basis of our analy- 
sis. Throughout the calculations, we adopt natural units with 
h = eo = c = 1. The electron mass is denoted by m. We 
make use of the Einstein summation convention unless stated 
otherwise, and we employ the length gauge for all wave func- 
tions and operators as we deal with off-resonant excitations, 
as it is done in most of the literature, and in textbooks on the 
subject (see e.g. |22]). The Mollow spectrum contains the in- 
coherent, inelastic part of the atomic fluorescence, i.e. the flu- 
orescence spectrum mediated by the many-photon processes 
whose intensity dominates over the elastic part in a strongly 
driven atom-laser system. In a purely quantum electrody- 
namic formalism, the description of many-photon processes 
would require perturbation theory in exceedingly high orders. 
However, as is well known, the description using dressed 
states [25] allows for a considerable simplification, as the for- 
mulas for the Mollow spectrum |24| follow rather naturally 
in terms of transitions among the dressed atomic states which 
incorporate the atom-laser interaction to all orders in pertur- 
bation theory. 

Before we now start with the discussion of the Mollow 
spectrum, a slight detour on questions related to gauge trans- 
formations of the laser-atom interaction is in order. The 
"length gauge" means that the laser-atom interaction is for- 
mulated in terms of the interaction —qEj_, ■ x, where q is 
the physical electron charge and Ei, is the (gauge-invariant, 
observable) laser field strength. Instead, in the "velocity 
gauge", the interaction is formulated as — g ■ p/m, where 
Al is a gauge-dependent, suitable vector potential for the 
laser field. In the velocity gauge, of course, one also has to 
add the term, but dipole interactions are mediated ex- 
clusively by the "Al • ^"-interaction. Due to gauge invari- 
ance, the two possible gauges are equivalent, provided that 
the gauge transformation of the wave function is properly 
taken into account E^- In order to avoid confusion, 
we stress here the absolute necessity of of considering the 
gauge transformation of the wave function in dynamical pro- 
cesses JllllllllllllllllllllllSiSIIIl. According 
to fl3ll2lll . the usual physical interpretation of a Schrodinger 
wave function is only conserved if the length gauge is used 
for the description of the atom-laser interaction. For dynam- 
ical processes, the velocity gauge leads to many more terms 
in intermediate steps of the calculation than the length-gauge 
formulation chosen here, due to the necessity of applying the 
gauge transformation to the wave function. Nevertheless, we 
would like to mention the possibility of an independent verifi- 
cation of our derivation, as presented here, in different gauges. 
In this case, the gauge transformation of the wave function 
should be applied already on the level of quantum mechanics 
(i.e., on the level of the Mollow spectrum as discussed in the 
current Section), not just to the quantum electrodynamic cor- 
rections discussed in the following Sections. This concludes 
our detour regarding gauge transformations. 

We recall from |24, 28, 29] that the incoherent resonance 
fluorescence spectrum of atoms driven by a monochromatic 
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coherent laser field may be expressed as 
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(2.2b) 
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of the driving laser field Ei^{t) = ^lel cos{uJi,t) with fre- 
quency cjl, macroscopic classical amplitude £l and polariza- 
tion El. q — —\q\ is the electron charge. A = o^l — <-^eg is 
the detuning of the laser field frequency from the atomic tran- 
sition frequency ojeg, and T is the spontaneous decay rate of 
the atomic transition. The excited and the ground state of the 
laser-driven transition are denoted by |e) and | g), respectively, 
and X is the position operator vector In secular approximation 
^ r, this expression simplifies to 



5inc(t^) 



r 
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Here, the separation of the Mollow spectrum into one central 
peak located atuj — oji^ and into two sidebands shifted by the 
generalized Rabi frequency JIr = VA^ + fj^ jj^^y easily be 
seen. The amplitudes and widths are given by 



Ao = 



A, 



02_h2A2 

2nl ' 



r± = r 



3J]2 _^ 2 A2 
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(2.5b) 
(2.5c) 
(2.5d) 



The approximate form does not represent the positions 
of the sideband peaks accurately in cases where F/i^R is not 



small. Indeed, the position uj± of the sideband peaks may be 
expanded in a series in powers of F/ whose first terms read 
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with y = A/0. For vanishing detuning A = 0, which implies 
Or = 0, Eq. (12. 6> specializes to 
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The correction terms move the sideband peaks closer to the 
central maximum. 

The above results in secular approximation may easily be 
interpreted with the help of the so-called dressed states, which 
are defined as the eigenstates of the interaction part of the 
Hamiltonian. Under the influence of the external driving field, 
the atomic states are no longer eigenstates of the Hamiltonian, 
but rather have to be combined with the driving laser field to 
give the new eigenstates. To show the precise composition of 
the dressed states, we use the quantum representations of the 
Rabi frequency 



= 2 5L\/n + 1 



(2 4) mixing angle 0„ defined by 

tan(26'„) = -0„/A , 
and the generalized Rabi frequency 



O 
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(2.8) 



(2.9) 



(2.10) 



rather than the corresponding classical entities. Here, n is the 
number of photons in the laser field mode, and the coupling 
constant for the interaction of the driving laser field with 
the main atomic transition is defined by 



9h 



-qigl^L ■ x\e)£j 



(7) 



(2.11) 



where = ^/lJiJ2V is the electric laser field per photon 
and V is the quantization volume. The matching of the electric 
field per photon with the corresponding classical macroscopic 
electric field £l is given by 



2 



If, 



(7) 



(2.12) 



Throughout this article, we will sometimes refer to the quan- 
tum description during the derivations, but use the classical 
entities in the final results. We may switch between the two 
descriptions as the driving laser field is assumed to be intense 
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in our analysis. The matching of the quantum and the clas- 
sical entities is possible with the help of the following list of 
replacements: 



A. Corrections to the Detuning 



1. Relativistic Corrections to the Resonance Frequency 



0n ^ 0, 

n + 1 « n . 



Using this notation, the dressed states are given by 



\{+,n)) ~cos6n \e,n) + sinOn \g,n + 1) 



(2.13a) 
(2.13b) 
(2.13c) 



(2.14a) 



-, n)) = — sin0„ |e, n) + cos6'„ \g,n+ 1). (2.14b) 



The well-known relativistic correction to the hydrogen en- 
ergy levels is given by 



H„ 



P 



nZa 



S{r) 



Za 



ctL. 



(3.1) 
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The effects are of the order of {Za)'^m, whereas the 
Schrodinger energy is of the order of (Za)^ m. The full ex- 
pression for the Dirac energy of a hydrogenic level with quan- 
tum numbers n, j, is |33, Eq. (2-87)] 



{Za)'^m {Za)*'m 
2v? ri^ 



1 



2j + l 8n 



, (3.2) 



Here, |i, n) (i £ {e, g}) are combined atom-field states where 
the atom is in state i with n photons in the driving field 
mode. The energies of these dressed states are given by 

E±^n = (n + V2)'^L + ± f2^''*/2, where the splitting 

between the two dressed states |(±, n)) is known as the AC- 
Stark-shift. The various spectral components then arise from 

transitions |(±,n)) |(±,n — 1)). The transitions H > + 

and > — yield the central Mollow component and the co- 
herent elastic peak, while the transitions H > — and > + 

yield the sidebands shifted to higher and lower frequencies, 
respectively. 



III. CALCULATION OF RELATIVISTIC AND RADIATIVE 
CORRECTIONS TO THE MOLLOW SPECTRUM 



In the following, we discuss corrections whose understand- 
ing is essential for the additional relativistic and radiative en- 
ergy shifts received by the dressed states. First, we evalu- 
ate corrections which may be incorporated in a redefinition 
of the detuning of the driving laser field to the atomic transi- 
tion frequency; in the second part, we complete the analysis 
by considering corrections which effectively modify the Rabi 
frequency. Throughout the analysis, we focus on the hydro- 
gen 1 S-2P transition as a promising candidate for a possible 
experiment. 

In Sec. mj we have employed a purely nonrelativistic the- 
ory. Both the resonance frequency as well as the transi- 
tion dipole moments are evaluated first for the nonrelativis- 
tic (Schrodinger) case. However, in order to resolve radia- 
tive effects, it is necessary to include the relativistic shifts 
of the transitions in a unified theory, and to analyze the 
fine-structure. The nonrelativistic expressions for the transi- 
tion dipole moments also change once we resolve the fine- 
structure levels, because the angular momentum algebra is 
augmented by the spin. 



where we neglect terms of order {Zaf . When evaluating the 
expectation values of H^ci on the dressed states ( I2.14l i in first- 
order perturbation theory, the following expression results. 



A 



Vf^2 + A2 



E, 



(3.3) 



Here 

eIT^ = (2^V2 l^reil 2Pv,) - {iSy, \H,,,\ 1%) , (3.4a) 
and 

= (2^3/, \Hroi\ 2P./,) - {iSi/, \Hroi\ ISy,) . (3.4b) 

The expression for Soj'^^'^j^ finds a natural interpretation as a 
first-order (in E^J ) correction to the quantity 



rol 



VfF + A2 



VfFTA2 



E. 



(3.5) 



We can thus formally define a "summed" relativistic shift of 
the Mollow sidebands as 

Au'^lf = ± (Jn^+(^A-El{l)' - + . 

(3.6) 

We recall that the detuning has been defined as A = cjl — 
in Sec. |lll If the fine-structure is included, the resonance 



■'eg 



frequency becomes j-dependent. The shift of the detuning as 
given in Eq. ( 13. 6> is thus equivalent to a modification of the 
resonance frequency according to 



^eg — ^eg ~t , 



A-K 



3l • 



(3.7a) 
(3.7b) 
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Thus, the "summed" shift of the detuning due to the relativis- 
tic correction E^.'^l , evaluated using the dressed-state basis, is 
equivalent to the shift of the detuning that would have been 
obtained if we had evaluated the detuning, right from the start, 
with a resonance frequency corrected by the relativistic ef- 
fects. The "summation" implied by Eq. ( 13.61 1 thus finds a nat- 
ural interpretation. 

Throughout the calculations, we will refer to shifts of 
the Mollow sidebands uj± due to first-order perturbations 
as Suj± [see e.g. Eq. (13. 3H . whereas summed expressions 
hke Eq. (13. 6> will be denoted 6lU±. 



of the detuning [in analogy to ( l3.7bH . 

With typical parameters (see Sec. llVt . the summed expres- 
sion Suj^.'^ yields results which significantly differ from the 

first-order expression Jw^"™^^ . The reason is that the bare 
Lamb shift is not small as compared to the detuning A, so 
that the higher-order terms of the series expansion are rele- 
vant. Nevertheless we use the summed formula Eq. ( I3.13t 
instead of Eq. ( 13. 8> . as it is the expected result in the sense 
that the Lamb shift is naturally interpreted as a modification 
of the transition frequency and a corresponding alteration of 
the detuning. 



2. Bare Lamb Shift 



3. Unified Expressions ft)r the Relativistic and Radiative Shifts 



In addition to the relativistic shifts, the positions of the side- 
bands have to be modified further if one desires a numeri- 
cal accuracy as required to appropriately model current high- 
precision spectroscopy experiments. In lEslEgil . the second- 
order radiative self-energy corrections due to the interaction 
of the combined system of atom and driving laser field with 
the surrounding non-laser-field vacuum modes was analyzed. 
Taking into account both interactions of the atom-field sys- 
tem with resonant and off-resonant intermediate states, in the 
limit A, f2 <C uj^g (i.e. under the replacements cjl ^ ^eg. 



— ljl ± f^R 0, and cjl 



"eg 



± JIr 'iujen) we obtain 



corrections to the energy of the dressed states which yield an 
additional shift of the position of the sidebands given by 



(Lamb) 



A 



tU) 
2^ bare ' 



(3.8) 



Here, the prefactor arises from the mixing coefficients cos 9 
and sin 9, and L\^^^.^ is the usual Lamb shift of the atomic bare 
state transition frequency which for the hydrogen 1 S-2P tran- 
sition is given by 



^hLc - ^2P, - Lis ■ 



(3.9) 



where j = 1/2, 3/2 is the total angular momentum quantum 
number of the excited state [for the definition of Lnij see also 
Eq. ( I3.16t below]. The Lamb shifts of the individual states are 
given by 134. .351 



Lis = 8172811(32) kHz, 
L2Pi^^ = -12835.99(8) kHz, 



•3/2 



12517.46(8) kHz. 



(3.10) 
(3.11) 
(3.12) 



The correction may be interpreted physically by defining the 
dressed summed Lamb shift 6uj±j as 




Vfi2 + A2 



r (Lamb) , 

(3.13) 

where to first order in ^^arc' recovers Eq. (13.8k Thus the 

(Laii 

±,j 



correction (Jw^^'J™'^'' effectively is a shift 



(3.14) 



Both the summed relativistic shift Eq. i3M as well as the 
summed Lamb shift Eq. ( I3.13> are effectively summarizing 
the corrections received by the detuning due to various shifts 
that go beyond the nonrelativistic treatment of the hydrogen 
(and Mollow) spectrum discussed in Sec. |lll These effects 
would also be observable in low-intensity scattering of (laser) 
light off atoms, and are automatically included in the ob- 
servable resonance frequency of the transition whose high- 
intensity behavior we are studying. The corrections can there- 
fore be included into the formalism if we replace the detuning 
A by the detuning Aoxp to the experimental transition fre- 
quency given by 



Acxp 



(3.15) 



where Wcxp is the experimentally observable transition fre- 
quency as it would be obtained from low-intensity scatter- 
ing I3al . In Sec. ini we have started from a nonrelativistic 
theory, and therefore the detuning A = ljl ^ t^eg was calcu- 
lated with regard to the inaccurate resonance frequency ujeg 
as it follows from the Schrodinger theory that fails, as is well 
known, to describe even the relativistic effects that lead to the 
fine-structure (let alone the Lamb shift). Thus, in practice, 
the bare Lamb shift modification to the detuning may be ac- 
counted for by replacing the resonance frequency ujeg as it 
would be obtained from a nonrelativistic theory, by an exper- 
imental value for the atomic transition frequency as found in 
low intensity scattering experiments 136,1 . 

The frequency Woxp may not be known well enough for any 
given transition to lead to a meaningful comparison between 
theory and a conceivable high-accuracy measurement of the 
Mollow spectrum. This is because we are sensitive, in the 
measurement of the Mollow spectrum, to tiny differences be- 
tween the laser frequency and the actual resonance frequency. 
It may therefore be useful to recall that for the Lamb shift Lni^ 
of a hydrogenic energy level (spectroscopic notation nlj), one 
may use the implicit definition (see e.g. 134. 31,] ). 

E{nlj) = nir [f{n,j) - 1] 

[/Kj) -1]' + ini, +^hfs, (3.16) 



2(m + mjv) 



where E is the energy level of the two-body-system and 
f{n,j) is the dimensionless Dirac energy, m is the electron 
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mass, rrir is the reduced mass of the system and rnjv is the 
nuclear mass. In very accurate experiments, one also has to 
include the hyperfine frequency shift E^s which depends on 
the quantum number toj? that includes the nuclear spin. Note, 
however, that the hyperfine structure does not contribute to the 
Lamb shift according to the definition Eq. ( I3.16> . 

The expression Eq. ( I3.16> can be used to make a theoret- 
ical prediction ojth for the transition frequency by forming 
the difference of this expression for the two states involved 
in the atomic transition whose high-intensity behavior we are 
studying. The detuning can then alternatively be evaluated 
as Ath = — Wth- Assuming wth = Woxp, one then has 
Apxn = Ath- For some recent data on Lamb shifts, we refer 

ton 



5. Off-Resonant Radiative Corrections Stimulated by the 
Atom-Laser Interaction 



For these corrections, we restrict the atom-field interaction 
to the laser mode, but take into account the off-resonant (OR) 
atomic levels \ j) (i.e. \j) ^ |e), \g)). The leading effect is the 
second-order perturbation 



(3.22) 



(±,n) 



E±,n — CHm + 'Hp) 



(±,n) 



4. Bloch-Siegert shifts 

The Mollow spectrum also receives corrections due to so- 
called counter-rotating interactions of the driving laser field 
with the atom fsS*]. These correspond to an excitation of the 
atom simultaneously with a creation of a laser photon or the 
vice-versa process. The first-order perturbation vanishes, and 
the second-order expression is given by |29j 



ae: 



(BS) 



= ± 



8 cos(26l) - (r^R/wL) [3 - 



1 

±- 



64 - 16 {Tln/iJLy 
2A2 + rj2 _^4AwL 



fi2 



A col 



(3.17) 



The correction Slu 



(BS) 



of the Mollow sidebands due to the 



Bloch-Siegert shift is thus given by 



6lu 



(BS) 



AE. 



(BS) _ 

A 



AE- 



(BS) 



(3.18) 
(3.19) 



Here, the parameter [B fi^] depends on the laser field intensity, 
which reflects the fact that the Bloch-Siegert shifts are stimu- 
lated processes. Assuming 51, A ^ ojl as a typical range of 
parameters, one has tJL ~ i^R and thus 



B 



1 



-0{A/ujln/ujl). 



(3.20) 



Here one should note that the energy shift Eq. ( I3.19> with B 
as in Eq. ( I3.20> vanishes for A = 0. This is consistent with 



the analysis in (291 l39il . where it was found that the Bloch- 
Siegert shift is suppressed by an additional power of for 
vanishing detuning. We also define a summed Bloch-Siegert 
shift in analogy with (13. 6> and (I3.13> as 



±(v/f^2 + (A -6172)2- 



A2) . 
(3.21) 

Effectively, the Bloch-Siegert shift may be accounted for by 
the replacement A ^ A — Bil^. This correction to the de- 
tuning is proportional to il^, i.e. proportional to the laser in- 
tensity. 



Here, we have defined the Schrodinger-picture Hamiltonian 
TYl describing the interaction of the atom with the driving 
laser field, the free energy of the nonresonant atomic states 
Ti-M, and the free energy of the electromagnetic fields Hf as 



Hl = -qx ■ El, 
Hm = 1-^') ^-^'l ' 



(3.23a) 
(3.23b) 

(3.23c) 



kA 



respectively, where CkA and a^^ are annihilation and creation 
operators for photons with wave vector k, frequency Wk and 
polarization A, and uoj (j g, e) are the energies of the non- 
resonant intermediate states. The sum in Eq. (I3.23c> extends 
over all possible vacuum field modes, and is the field op- 
erator for the laser mode. 



Et 



(3.24) 



Here, is the polarization vector for the laser mode. As the 
laser mode is highly populated with an occupation number 
n ^ 1, both the field annihilation and creation operators in 
TCl contribute [see also Eq. ( I2.l2l i1. The resulting expression 
for the energy shift of the dressed state |(+, n)) is given by 



n + I n 



-ujj + El -ujj + E2 
n + 1 n 



|5,,f sin^^ 



-W,- + i?3 — W,- + £'4 



(3.25) 
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with 



3 



UJg + 2 WR, - 

w„ + 2 wr , 



Ea — UJn 



- WR, , 



Wr. 



(3.26a) 
(3.26b) 
(3.26c) 
(3.26d) 



For the dressed state | (— , n)), we have 



(OR) 



E 



\9ej\^ sin^i 



n 



1 



, + E, 
1 



-t^j + E3 



—UJj + i?4 



(3.27) 



In calculating these expressions, we may carry out the semi- 
classical approximation n + 1 « n and replace the above Rabi 
frequency by its semiclassical counterpart. The coupling gij 
is defined by 



(3.28) 



and is of the same order of of magnitude as We therefore 
obtain as the second-order shift due to the off -resonant energy 
levels. 




as: 



(OR) 

±,Tl 



AK 



(OR) 



±V cos(26l) 



where A is the detuning and 



A 



V 



\9e]\ 



+ 



[V^f], (3.29) 



1 



E; 



El —UJj 
1 

-U^, + Ea 



Eo 



(3.30) 



depends again on the laser intensity. The energies Ei (i = 
1. . . . , 4) are defined in Eq. (I3.26> . Therefore this additional 
shift is a stimulated radiative correction in the same sense as 
the corrections discussed in the previous section fy3\. 

To further evaluate the parameter T), it is important to note 
that the virtual states are coupled to the initial states by the 
driving laser field rather than by the vacuum as for example 



in Lamb shift calculations. Thus the polarization of the cou- 
pling field mode is fixed. For the off-resonant corrections, it is 
sufficient to evaluate the relevant matrix elements in the non- 
relativistic approximation. So, if we assume the atoms to be in 
the 1 S ground state initially (to = 0, of course), then we have 
a situation in which for a given polarization of the laser field 
not all magnetic sublevels of the 2P states are coupled. In the 
following, we assume the driving laser field to be z-polarized 
[cl = (0, 0, 1)], so that only the m — sublevels of the IS* 
and the 2P ground and excited state are occupied. Then the 
parameter V may be rewritten as [see Eq. ( I3.26l l1 



M{9\^\e)\' 

in terms of the two matrix elements 



MgiO - (^g\zG"iOz\g) 



MeiC) 



■■G (C)^ 



where the propagator is given by 



G(C) 



1 



(3.31) 

(3.32) 
(3.33) 

(3.34) 



and where the double prime means that both resonant states 
\g) — jlS*, TO = 0) and |e) = |2P, ?7i = 0) are exc/Mt/ec/ from 
the sum over intermediate states in the Green function. The 
evaluation of Afe(C) requires special care, as there are both 
S (angular quantum number I = 0) and D (I = 2) states as 
intermediate states. Due to the fixed polarization of the cou- 
pling field, the angular parts of these two contributions have 
different proportionality factors relative to the angular parts of 
the "standard" matrix element 



(3.35) 



and thus have to be calculated separately. The above standard 
matrix element may however be recovered from the matrix el- 
ements for definite initial state and coupling field polarization 
by averaging appropriately. This is discussed in AppendixlEl 
In the propagator G(C), the energy is parameterized by the 
(in general complex) parameters 



c ^ m 



[Zafm 
2n2 i2 



Za 
n 



TO 



(3.36a) 



(3.36b) 



Usually, one has ( — E—lu where E is the bound state energy, 
and we may write 



at) = E- Lo{t) 



{Za) 



2n2 



uj{t) 



(337) 
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The parameters lo and t are related by the equations 



2v? t2 

1 



m(Za) 



-1/2 



(3.38a) 



(3.38b) 



where n, / are the principal and the angular momentum quan- 
tum number of the quantum state for which the relevant matrix 
elements are to be evaluated. In the following, we will also use 
common spectroscopic notation for the level characterized by 
n and I, i.e. for example t2p{Lo) = t2i{Ld). For the energies 
i^i (i € {1, . . . , 4}) and n = 1, 2 we thus obtain 

(3.39) 



The above matrix elements Eqs. (l3.32ll3J3l without the dou- 
ble primes, i.e. including resonant intermediate states, may 
then be expressed in terms of the "standard" hypergeometric 
function 1 34 .4(1 4 1.1 



t) = 2F1 1, -nt, 1 - 7it, 



as (see AppendixlEl 



Mg{0^{9\zG{C)z\9)^ni4 



1-r 

i + t. 



3(t- + 



256 



3(^-1)5 [t + lf 
Me(C) ={e\z G{C) z I e) = m a| 
2i4<ii(23t2-7) 



(3.41a) 



15(t-l)7(t+l) 



(3.41b) 



where 



Xg{t) ^ 38t^ + 26t^ + - m* - 12t^ 

+ I2t^ + 3t - 3 , 
Xe{t) = 6739t^" - 1702t^ - 231t® - 1420t^ 



262r 



1944t^ - 402i* - 1140t^ 



435i2 + 270t - 135 . 



(3.43) 



Here, the Bohr radius scaled by the nuclear charge number Z, 
in our units, is 



1 



Zam 



(3.44) 



where a is the fine-structure constant and m is the electron 
mass. 

The corresponding matrix elements without resonant inter- 
mediate states may then be obtained by subtracting the respec- 



tive contributions of the resonant intermediate states 

M,(C) = {g\zG"{0^\9) 



MAO '^''^'^^ 



E2P — c 



Me(C) 



zG {Qz 



MeiC) 



l(g|^|e)r 



(3.45) 



(3.46) 



We assume here that the Rabi frequency is not excessively 
large, which implies that it is small as compared to the op- 
tical transition frequency (i.e., ^ wr). For a meaningful 
measurement of the Mollow spectrum, it is necessary, further- 
more, to tune the laser close to the atomic resonance (which 
implies A ^ ojr)- In this case we may carry out the following 
approximations [c.f. Eq. ( I3.26H 



El 

E2 

E3 
Ei 



2wR 



t2p{Ei) 
t2p[E2) 

hsiEs) - 



■ 1/2 , (3.47a) 

•V^Ol, (3.47b) 

2/V7, (3.47c) 

2 . (3.47d) 



(3.40) With these parameters, V evaluates to 



V 



1 



(Za)'^m 



[6.2148(8) - 0.23532(2)1] 



— [2.3305(3) - 0.088245(6) i] 



(3.48) 



The uncertainties are mainly due to the approximations car- 
ried out in Eq. ( I3.47> with respect to the energies Ei (i — 
1, . . . , 4) originally defined in Eq. ( I3.26t . This is possible be- 
cause the off -resonant stimulated radiative correction amounts 
to a modification of the detuning which is of order ft^/uij^, <C 
[see also Eq. ( I3.49> below]. Therefore we may carry out 
the approximation Eq. ( I3.47t . i.e. neglect the further correc- 
tions of order fl^/uj^ ^ fi^/t^R and fi^A/w^ <C Q'^/ujb., 
which are beyond the scope of the current analysis. The real 
part Pr = Re (P) gives rise to a shift of the position of the 
Mollow sideband, while the imaginary part = Im (P) de- 
scribes the ionization into the continuum. This means that the 
(3.42) imaginary part of the energy shifts AE^^"^ received by the 
two dressed states |(±, n)) must be negative, which is equiv- 
alent to a negative sign for the imaginary part Pj. 

In the numerical analysis in Sec. IIVI it is shown that for 
typical parameters the system is sufficiently far from the ion- 
ization threshold |42] so that the ionization does not restrict 
the applicability of our scheme. The real part yields a correc- 
tion to the detuning given by A ^ A — PrJ^^, according to 
the summation [cf. Eq. ( I3.29H 



6lu 



(OR) 



± ( y^n^ + (A - Pr1]2)2 _ ,/q2 ^2 



with 



Vf^2 + A^ 



(3.49) 



[pR n^] . (3.50) 
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As pointed out below in section irVl (see also Tab.|IJ, the mag- 
nitude of the off -resonant correction is small as compared to 
the detuning for typical parameters so that there is no numer- 
ically significant difference between the first-order correction 
Eq. ( I3.29> and the summed form given in Eq. ( 13.491 1. 

The first-order imaginary contribution to the Mollow side- 
bands is 



(3.51) 



This effect broadens the sidebands slightly, but its contribution 
is so small for typical parameters (see Sec. lIVI below). that it 
may be ignored on the current level of accuracy of the theo- 
retical predictions. It is interesting to note that the stimulated 
off -resonant correction is small as compared to the relatively 
large effect mediated by the bare Lamb shift of the transition, 
which is discussed in Sec. IIII A"2l 



B. Corrections to the Rabi Frequency 

In this section, we discuss corrections to the incoherent 
fluorescence spectrum due to modifications of the Rabi fre- 
quency. In particular, we consider corrections to the transition 
dipole moment and to the spontaneous transition rate of the 
atomic transition, due to coupling of the driving laser field to 
resonant and nonresonant atomic transitions, and corrections 
to the secular approximation. 



1. Relativistic Corrections to the Transition Dipole Moment 

In this section we discuss relativistic corrections to the flu- 
orescence spectrum up to relative order (Za)^. The correc- 
tions amount to a modification of the atomic transition fre- 
quency and of the transition dipole moments. The relativis- 
tic expressions for the state energies and the transition dipole 
moments depend on the total angular quantum numbers j of 
the involved states, which is the vectorial sum of the electron 
orbital angular momentum / and its spin s. Therefore we spec- 
ify the total angular momentum quantum number and thus the 
spin state of the atomic system in order to fix a specific ex- 
perimental setup (see also Appendix We further assume 
the atom to be driven by a pure dipole laser field linearly po- 
larized in z-direction, such that the laser field only couples 
states with equal magnetic quantum number. The situation of 
a pure dipole field has recently been studied in a related con- 
text in 143]. In the numerical analysis in Sec. II VI we consider 
a standing wave laser field configuration where the atom is 
at a point of maximum electric field intensity of the standing 
wave. As then the magnetic field component of the driving 
laser may be neglected, it is not considered in the following 
analysis (corrections due to the variation of the electric field 
about its maximum are treated in Sec. IIII B"2t . 

The relativistic corrections to the energies of the atomic 
states and thus to the atomic transition frequencies effec- 
tively modify the detuning A and may be accounted for by 



choosing an experimental value for the atomic transition fre- 
quency as found in low-intensity scattering experiments (see 
Sec. IIII A"3l . The corrections to the transition dipole moments 
may be evaluated with the help of the relativistic wavefunc- 
tions of the hydrogen atom as given in |44, 45J. We denote 
the absolute relativistic correction to the nonrelativistic ma- 
trix element 



(lS'i/,,m = ±V2| z |2Pj,m = ±V2)r 



(3.52) 



by ddj , where j — 1/2, 3/2 is the total angular momentum 
and TO — ±1/2 is the magnetic quantum number. In the fol- 
lowing, we will omit the "to =" from the second parameter of 
the atomic state vectors. Then the relativistic matrix element 
(with subindex "R") gives rise to a relative (ZQ;)^-correction 

SS^^ /dj with respect to the nonrelativistic (NR) expression 
which is given by 



(15v„±V2|z|2P„±V2)r 
(1%,±V2U|2F,,±V2)nr 



= 1 



5d\ 



(R) 



(3.53) 



where we ignore higher-order relativistic terms of order 
(Za)™ with TO > 4. The corresponding matrix elements 
where the "initial" and the "final" state have different mag- 
netic quantum numbers vanish identically as the driving laser 
field is assumed to be polarized in the z-direction. Evaluating 
the relative corrections, one obtains 



""1/2 



di 



Sd. 



(R) 

i/2 



da 



:/2 



[Za] 



{Za) 



13 
32 

31 
96 



In2-ln3 



ln2- 



■ln3 



(3.54a) 



(3.54b) 



The Rabi frequency and the transition dipole moment depend 
linearly on each other Therefore the relative correction to 
the Rabi frequency is identical to the relative correction to the 
transition dipole moment: 



5Vl 



(R) 



6d 



(R) 



(3.55) 



Here, 5Vlj is the absolute correction to the Rabi frequency 



due to relativistic modifications to the transition dipole mo- 
ment. The spin-dependent shift (5tU^j of the position of the 
Mollow sidebands due to the relativistic corrections of the ma- 
trix element is thus given by: 



where 



and 



,,,_C + Ora| . (3.58) 



(3.56) 



: £i (3.57) 
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With these definitions, the summed relativistic correction to 
the dipole moment effectively corresponds to a replacement 

£j). 



This naive expansion gives the wrong prefactor and is not ap- 
plicable here. 

The term in Eq. ( I3.61> entails a spin-independent correction 
to the transition dipole moment. At resonance, one has 



2. Field-Configuration Dependent Correction to the Rabi 
Frequency 

It is well-known that the magnetic component of plane- 
wave electromagnetic wave influences the transition current 
at relative order {ZaY (see e.g. |34, 0, S). The radia- 
tion pressure due to the magnetic field could move the atom. 
Therefore, we propose a standing-wave field configuration, 
where the atom is placed at an anti-node of the standing-wave 
electric field. In this setup, the influence of the magnetic field 
can be neglected to a very good approximation. 

The analysis of the previous Sec. IHIBTl is valid up to the 
order discussed (relative order (Za)^) only for a pure dipole 
field which additionally has to be constant in any direction 
perpendicular to the polarization. However, for a standing 
wave configuration, the z-polarized electric field of the laser is 
not constant in the propagation (a;-)direction. This leads to a 
further correction, which gives rise to a field-configuration de- 
pendent shift of the Rabi frequency. In the following, this shift 
of relative order {ZaY is analyzed for the setup described 
above where the atom is at the maximum of the standing-wave 
electric field. 

We start from the long-wavelength quantum electrody- 
namic (LWQED) interaction Hamiltonian |47|. The only rel- 
evant terms (in the context of our analysis) of the interaction 
part of this Hamiltonian are 



-qx-E- j-x' E\ 



6 



■ x' x^ E\ 



(3.59) 



Here, the a;* denotes the ith component of the position opera- 
tor vector X, and E^j is the partial derivative with respect to x^ 
of the ith component of the electric field vector The electric 
field of the standing wave is given by 



E{t, x) — iz Ssw cos(ijji) cos{kx) . 



(3.60) 



The term containing the first derivative of the electric field in 
Eq. ( I3.59> vanishes, and the last term gives 



^ i j k T?i 



' X X X Ej a V, ^ z X El ^ ^ 

6 6 ' 



^zk^ x^ £sw cos(tjt) . (3.61) 



This result has to be distinguished from a simple expansion of 
the electric field around the maximum at a; = 0, which yields 

— qr ■ E{t,x) — —qz £sw cos(cjt) cos(fca:) 



-qz £sw cos(wi) I 1 



{kx)"' 



0{x^) 



(3.62) 



k = - {Zafm. 



(3.63) 



The relative correction due to the additional contribution 
Eq. ( I3.61> to the interaction Hamiltonian is therefore 



fc2 (lS',m = 0|zx2|2P,m = 0) 
Y (15,TO = 0|z|2P,m = 0) 



(3.64) 



Analogous to Eq. ( I3.53> . this modification of the transition 
dipole moment, for a 15-2^ transition, gives rise to a correc- 
tion to the Rabi frequency given by 



1 

16 



{Zafn. 



(3.65) 



The summed shift 5lJj\ ' of the Mollow sidebands due to this 



modification of the Rabi frequency can be expressed as 
where 



(3.66) 



vn2 + A2 



and 



16 



(Za) 



(3.67) 



(3.68) 



With these definitions, the summed relativistic correction to 
the dipole moment effectively corresponds to a replacement 

It is important to note that the long-wave QED correction 
to the interaction Hamiltonian Eq. ( I3.59t does not couple any 
unwanted magnetic quantum numbers to the laser-driven dou- 
blet. 



3. Higher-Order Corrections ( in O, and A) to the Self-Energy of 
Dressed States and Corresponding Correction to the Rabi 
Frequency 

The Lamb shift of dressed states is different from the Lamb 
shift of atomic bare states, as already discussed in SecU In 
this section, we extend the analysis of Sec. lIII A'2l to the next- 
higher order For this, we keep the terms linear in and A in 
evaluating the energy shifts of the dressed states. As explained 
in detail in [28, 29J, we thereby obtain a further correction 
to the position of the fluorescence sidebands which may be 
expressed as 



(3.69) 



11 



where 



C 



ln[(Za) 



(3.70) 



is a dimensionless constant. For the hydrogen 1 S-2P transi- 
tion, the leading logarithmic term is independent of the spin 
and given by 



C = —aiZaf \n\iZa)-^] 
An 



(3.71) 



This corrections may be interpreted physically with a summa- 
tion as used for the bare Lamb shift correction: 

sJf^ = ± (y/riP~(r~CyTA^ - VQ^ + aA , (3.72) 



(C) 

6u! j. because of the smallness of the correc- 

(C) 



with Suj 

tion. Thus the additional shift Slu)^' may be interpreted as 
radiative correction fl 17(1 — C) of the Rabi frequency. 



4. Leading Logarithmic Radiative Corrections to the Transition 
Dipole Moment (Vertex Corrections) 

In "normal" bound-state quantum electrodynamics, vertex 
corrections are evaluated with respect to the interaction of the 
electron with the binding field of the atomic nucleus. In an 
effective treatment, and in leading logarithmic approximation, 
the effect of the self-energy may be accounted for by making 
use of an effective Lamb-shift potential |48] 



Aliamb(r) ^ -a{Za)\n[{Zay 



m 



which modifies the Coulomb interaction according to 



Za Za , . 

> +AVLamb(r-). 

r r 



(3.73) 



(3.74) 



Note that the potential j3.73t is really the consequence of a 
self-energy ("vertex") correction, not that of vacuum polariza- 
tion. In many cases, vacuum polarization corrections may also 
be accounted for by employing an effective potential, but the 
corresponding potential lacks the large logarithm ln[(Za)^^]. 
An accurate treatment of self-energy corrections requires the 
consideration of many more terms than the crude approxima- 
tion ( I3.73> . 

Here, we evaluate the leading vertex corrections to the in- 
teraction of the bound electron with the driving laser field, 
and so we have to consider both the Coulomb as well as the 
laser field. Nonrelativistically, the atom-laser interaction is 
given by the matrix element of the usual interaction Hamil- 
tonian, which reads (length gauge, 15* <^ 2P transition) 
-qEL{lS\z\2P) = -qEi^d where is the field strength 
of the (strong) laser field and the dipole moment is 



{2P\z\lS) 



-V2^— 

3^ m{Za) 



(3.75) 



Vertex corrections lead to modifications of the dipole moment 
given by d ^ d + Sd, where the vertex correction 6d is 
considered below, and the radiative correction (in the length 
gauge) to the laser-atom interaction is effectively a replace- 
ment —qEi^ d —qEi^ [d + &d). The large intensity of the 
driving laser field is accounted for in this formalism because 
the electric laser field strength i?L multiplies both the dipole 
moment matrix element of the interaction Hamiltonian d and 
the radiative correction 5d. 

Laser photons as well as the spontaneously emitted photons 
in the radiative decay of excited states are real rather than vir- 
tual. Consequently, the radiative corrections to the laser-atom 
interaction on the one hand and to the radiative decay rate on 
the other hand are related to each other In the length gauge, 
the leading-order expression for the spontaneous emission de- 
cay rate is 



r 



-a {Za)^ m . 



(3.76) 



In order to obtain gauge-invariant results for the quantum 
electrodynamic corrections to F (while working in the length 
gauge), it is necessary to consider both radiative vertex correc- 
tions to the dipole moment d and corrections to the (bare) tran- 
sition frequency (energy difference) E (see Refs. |48, 49, 50|). 
In our treatment, the vertex corrections to the dipole moment 
are given in Eqs. (I3.81> . (|^6^ and ( l3.86bL whereas the ver- 
tex corrections to the transition frequency enter into the radia- 
tive correction to the detuning in Eq. (I4.10a> . 

In general, the vertex corrections to the laser-atom inter- 
action enter at the relative order of 0{a{ZaY 1ii[(Zq;)^^]). 
One may wonder why the corrections do not enter at the rel- 
ative order 0{a). The reason is that in interactions with real 
photons (the square of the four-momentum being q^ = 0), the 
otherwise dominant correction due to the Fi Dirac form factor 
vanishes, and the remaining terms are then of higher order in 
the Za-expansion. 

We now analyze the shift of the resonance frequency and 
the shift of the transition dipole moment induced by the Lamb- 
shift potential ( I3.73> . The transition energy 

E = E2P- Eis = I [Zafm (3.77) 
o 

is shifted by VLamb according to 

^ijdog) = -(l5|AFLa,„b(r-)|15) 

= ^—a{Zaf\n[{Za)-^], (3.78) 

because the matrix element of the 2P state vanishes. This 
yields a relative shift of 



(5£;(i°8) 32 



E 



^a{Zaf ln[(Za)"^] 



(3.79) 



The modification of the matrix element due to the corrections 
to the 1 S wave function amounts to 



= ( 2P 



1 



E-H 



AVr 



Lamb 



(r) 



15 



%/2^M^ln[(Za)-2]f481n^ + 13l), (3.80) 
■i'rmr \ 6 J 
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where the prime denotes the reduced Green function. Thus 
the logarithmic relative correction is L49il 



— — = ^ \n[{Za) 2] 



This leads to a correction to the Rabi frequency analogous to 
Eq. (l335l . 



(3.82) 



where we may ignore the spin, in contrast to ( I3.54> . 

The interpretation of this shift is analogous to the relativis- 
tic corrections to the dipole matrix element in Sec. 1111 BTI The 
MoUow sidebands are shifted by the frequency 



5uj 



± 



± 



where 



and 



A 



±- 



-4, 



O 



A2 
(3.83) 

(3.84) 
(3.85) 



is a dimensionless constant. Then these corrections may be 
accounted for by the replacement ft — > il{l+A). One should 
note that at this parametric order the above results hold for 
both states 2Pi/^ and IPa/^ I49. „,5a]. 



5. Nonlogarithmic Vertex and Vacuum Polarization Corrections 
to the Transition Dipole Moment 

The concurrence of the radiative shifts to the transition 
dipole matrix elements for the 15 — 2P1/3 an d the 15 — 2P3/2 
transitions, which was found in Sec. 1111 B"4l for the effects of 
relative order a{Zaf ln[(Za)^2]^ jg lifted on taking into ac- 
count corrections of order a{Zaf' (no logarithms) due to the 
self-energy. Vacuum polarization corrections to the transi- 
tion dipole matrix elements also enter at the relative order of 
a[Zaf. 

The sum of the nonlogarithmic ("nlog") vertex and vacuum 
polarization corrections to the transition dipole moments of 
order a (Za)"^ are given by isoll 



5d\ 



(nlog) 



a{Za) 



di/2 

for the 2P1/2 state and by 

d3/2 



9.2(1. 



(3.86a) 



a{Za) 



9.3(1.9) 



(3.86b) 



for the 2P3/2 state (sd]. In between Eqs. (70) and (71) of fs^, 
it is stated that currently, there is an internal discrepancy be- 
tween the numerically obtained values for the radiative cor- 
rection at Z = 5, 10, ... on the one hand and analytic results 
for the first terms of the Za-expansion (logarithm+constant) 
on the other hand. This discrepancy is of the order of 10 % 
of the total constant term of order a{Za)^, and this limits the 
current status of the theory. Here, we employ an even more 
conservative error estimate and assign a 20 % uncertainty to 
both of the numerical values in Eqs. j3.86a> and j3.86bt . The 
corresponding correction to the Rabi frequency is again anal- 
ogous to Eq. ( I3.55> . 



6n 



(nlog) 



6d 



(nlog) 



(3.87) 



One obtains the following total correction SlJ^^^^"^ to the po- 
sition of the sidebands due to logarithmic and nonlogarithmic 
correction to the transition dipole matrix element (TDM): 



60J 



(TDM) _ 



± 



where 



6uj 



(TDM) 



r!2 



and 



5d^'°s) Sd' 



(nlog) 



d, 



0(a{Za) 



+ A2 
(3.88) 

(3.89) 
) . (3.90) 



Aj is a dimensionless constant and j = 1/2,3/2. Thus the 
corrections amount to a modification of the Rabi frequency 
given by fl -> $1(1 + Aj). 

We expect a similar nonlogarithmic correction of order 
a{Za)'^ to supplement the C-term discussed in Sec. 1111 B 31 
The logarithmic C-term was found to be of the order of 
a(Za)2 ln[(Za)^2] xhe evaluation of the corresponding 
nonlogarithmic term is however beyond the scope of this 
work. Here, we only present a (conservative) estimate of the 
expected correction. To this end we observe that the general 
form of the radiative corrections under discussion is 



a{Za)'^ ln[{Za) ^] ■ Ci + a{Za)'^ ■ C2 



(3.91) 



where ci and C2 are dimensionless constants. Based on ex- 
perience with similar corrections (both for self-energy effects 
and well as radiative corrections to decay rates [37, 50]), we 
assume the following relation with the corresponding uncer- 
tainty for the unknown parameter C2'. 



C2 



-2± 2) • ci 



(3.92) 



For example, we verify the validity of this estimate for the 
particular contribution to the nonlogarithmic part discussed 
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above, i.e. for the terms that contribute to Aj . For these terms, 
the estimate (—2 ± 2) • ci evaluates to —2.56 ± 2.56 for the 
constant term. In comparison, the values obtained above for 
C2 by a direct numerical analysis are —2.92 for the correc- 
tion Ai/2, and —2.97 for the correction .43/2, which agrees to 
the estimate. For the C-term correction, we thus obtain from 
Eqs. (m\ and (i32\ 



a{Zaf ( ^ ln[(Za)-2i 



47r 



(2 ±2) 



(3.93) 



as the combination of the logarithmic and the estimated non- 
logarithmic correction. We expect the nonlogarithmic correc- 
tion to be spin-dependent in analogy to ( I3.86> . As for the cor- 
rection to the matrix element, this shifts the Mollow sidebands 
by [see also Eq. ( I3.72H 



(3.94) 



such that the correction may be applied by the replacement 



6. Corrections to the Secular Approximation 

In this section, we transform the correction terms to the sec- 
ular approximation in Eq. \2.b\ such that they may be inte- 
grated into our correction scheme. Eq. ( 12. 6> may be rewritten 
as 



f]2 



2 /T^\ 4 



5(1 + 2/2) \n) ^""yn 



(3.95) 
(3.96) 



For A <C il, one may expand the leading contribution to S to 



give 



O 



(3.97) 



With this definition, the corrections to the secular approxima- 
tion may be accounted for in the final result with the replace- 
ment £7 — > (1 — 5), which results in a summed shift of 



6u! 



,(S) 



±(v/r!2(l-5)2 + A2- ^ 



(3.98) 



IV. NUMERICAL DATA FOR THE HYDROGEN 1S-2P, 
TRANSITIONS (j = 1/2, Vz) 

In the previous section, we have discussed both corrections 
to the detuning and to the Rabi frequency, which give rise to 
a modification of the positions of the sidebands in the Mollow 
spectrum. We start the evaluation of correction terms here 



from a point where we assume that all relativistic corrections 
to the transition frequency, as well as hyperfine-structure ef- 
fects, have already been included in the bare transition fre- 
quency. This frequency corresponds to the prediction obtained 
using Eq. ( I3.16> by setting L„;^ explicitly equal to zero. We 
therefore redefine the detuning A to be 



where [cf. Eq. ( I3.16l l1 

= EC^Pj) - EilS) 

with 

E{nlj) = nir [f{n,j) - 1] 



(4.1) 



(4.2) 



2(m + m^v) 



hfs • 



(4.3) 



Of course, the full theoretical prediction E{nlj) is obtained 
as the sum of E{nlj) and Lni^- The modification of the Rabi 
frequency due to the bare Lamb shift 



(4.4) 



as discussed in Sec. lIII A"2l finds a natural interpretation as a 
contribution to the Lamh shift of the dressed states. The vari- 
ous first-order correction terms to the lowest-order prediction 
for the generalized Mollow-sideband displacement 



±Vf^^+A2, 
starting from (14. 1> . may be summarized as follows: 



(4.5) 



A 



• T U) 
^2 \ bare 



±- 



" Vf72 + A2 



where 



-(1/2) ^. 

''bare 
''bare 

B 



8185.647(34) • 10^ Rz, 
8160.294(34) • 10^ Uz . 



A 

A1/2 
A3/2 



1 

4cjR, ' 

X'=X>R, + il?i 

_ 1 

a(Za)2 



2.3305(3) - 0.088245(6)1] 



=A 



a{Zaf 



9.2(1. 



--A 9.3(1.9), 



(4.6) 

(4.7a) 
(4.7b) 

(4.7c) 

(4.7d) 
(4.7e) 

(4.7f) 
(4.7g) 
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C,^a{Zaf—{ln[{Zar^]~{2±2)) , 



S = 



-(- 



(4.7h) 
(4.7i) 

(4.7j) 
(4.7k) 
(4.71) 



We have found that the bare Lamb shift L^^^^, the Bloch- 
Siegert shift (B), and the off -resonant self-energy corrections 
(I?r) give rise to a modification of the detuning A in the ex- 
pression for the Mollow spectrum. As the latter two effects 
are intensity-dependent, also their correction to the detuning 
depends on the intensity of the incident laser field. The modi- 
fications to the transition dipole moment (Aj , £j , T) as well as 
the higher-order resonant self-energy shifts (Cj) and the cor- 
rection to the secular approximation {S) may be interpreted as 
radiative corrections to the Rabi frequency VI. The interpreta- 
tions as a modification of the detuning and the Rabi frequency 
may best be seen by using a summation of Eq. (14. 6> . which 
we have shown to be valid up to first order in the parameters 
in Eq. dO : 



with 



(4.8) 



(4.9) 



Here, the Rabi frequency and detuning are supplemented by 
the discussed relativistic and radiative corrections; these are 
given by 



aU) _ tU) 

rad bare 

0) 



rad 



■'4,') ^ 1 ^ 1 



T ~S. 



(4.10a) 
(4.10b) 



This summation implied by ( 14. 8> is motivated by Eqs. ( 13. 6> . 
( |3JH , dm , ( IT49t , ( l336t , ( IT72t , dl88t , dl66l , dlWl , and 
( I3.98> . In (l4.10bL the symbol f]J;^|[ indicates a relative modi- 
fication of the Rabi frequency, i.e. a dimensionless quantity. 

We have thus "summed" all the radiative corrections as 
effective corrections to the Rabi frequency and the detun- 
ing. Of course, the mixing angle as defined in Eqs. ( 12. 9> 
and (|^2D is changed by the radiative corrections. Indeed, 
one may evaluate the corrected 9 by employing the relation 
tan(2e') = -ficorr/Acorr whcrc Acorr = A - a|,^^j and 
ficorr = • ^1 + ^^rad) rclativistically and radiatively 

corrected Rabi frequency and the detuning, respectively [see 
Eqs. (^^bj and ( l4.10aH . Because all relativistic and radia- 
tive corrections find a natural interpretation as corrections to 
the Rabi frequency and the detuning, the corrected dressed 



states have the same structure as Eqs. ( I2.14a> and (I2.14b> . but 
with relativistic and radiative wave function corrections a cor- 
rected mixing angle. 

All relativistic and radiative corrections to the Rabi fre- 
quency and the detuning have been evaluated here using the 
unperturbed mixing angle 6. The "corrections to the correc- 
tions" (sic!) due to an evaluation of modifications to the Rabi 
frequency and the detuning in terms of the corrected mixing 
angle are of higher order than the terms relevant for the dis- 
cussion in the current paper and may be neglected on the level 
of approximation employed in the current investigation. 

In the following numerical analysis, we assume the atom 
to be located at an anti-node of a laser field in standing-wave 
configuration. The atom is thus driven by two counterprop- 
agating laser beams, whereas in the definitions of the elec- 
tric field and the Rabi frequency in Sec.|n|and especially in 
the matching of the classical macroscopic field with the cor- 
responding quantum counterpart in Eq. \2.\2\ a single-mode 
running-wave field was considered. However all results of 
Sec. Hn] also apply to a standing wave field configuration if 
the macroscopic electric field strength £l is taken to be the 
total field strength of both counterpropagating field modes at 
the position of the atom. In order to avoid confusions, from 
now on we thus denote the total electric field strength of the 
standing wave as f sw- 
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-/2 



2P, 



-/2 



For the \2Pi/.2,m = ±1/2) state as upper state, the decay 
constant is given by Ti/^ = 99.70942(1) • 10^ Hz |50]. In 
order to account for the dependence on the laser field intensity, 
we introduce the parameter ^,1/2 = |ri|/ri/2. Then for hi/.^ = 
1000 and A 50 Ty^, one has wr > O > Ti/a, A. 
Therefore, the relative corrections to the detuning and the Rabi 
frequency in Eqs. ( 14.1 Oat and ( I4.10bt become 



rad 



-8.175249(33) • lO^Hz, (4.11a) 
-19.78(56) • 10"^ (4.11b) 



Here, the parameter /11/2 may be expressed in terms of the 
electric field strength £l as 



hi/,_ = 346.783 • 10"^ |£sw [V/m]| 



(4.12) 



Of course, [V/m] in this case means that the peak electric 
field strength of the laser is assumed to be measured in Volts 
per meter. 

The absolute ionization rate I into the continuum due to 
the driving laser field is given by X = ViVP'. In an experi- 
ment, this ionization rate has to be much smaller than the Rabi 
frequency, such that on average the atom undergoes many flu- 
orescence cycles before it is ionized. Thus we define the ratio 



In = 



(4.13) 
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Shift 


ISi/a ^ 2Pi/2 [kHz] 


IS 1/2 <-* 2P3/2 [kHz] 


^_(Lamb) 


741599(4) 


738281(4) 




-50.30(5) 


-50.30(5) 




-468.51(6) 


-468.51(6) 




-1842.1(1) 


-1937.7(1) 




-331.44(2) 


-331.44(2) 




-121(31) 


-121(31) 




374(25) 


372(26) 




-49.8(2) 


-49.8(2) 



TABLE I: Summary of all individual energy shifts due to the vari- 
ous discussed corrections. All numbers are obtained for hj = 1000 
and A = 50 ■ Tj with j = 1/2 (j = 3/2) for the left (right) column 
using summation formulas such as Eq. <4.8> . Here, (5Zj*|^^™'^' is the 
correction to the high-frequency MoUow sideband position related 
to the bare Lamb shift (c.f. Sec. IIII A2> . the symbol (BS) denotes 
Bloch-Siegert shifts (Sec. IIIIA4> . and the (OR)-shifts are due to 
off-resonant excitations (Sec. IIII A 5> . These shifts all may be inter- 
preted as arising from a modified detuning A and are discussed in 
Sec. lIII a1 The other five corrections are due to a modified Rabi fre- 
quency f Sec. IIII B1 . In particular, SZj'^^ is discussed in Sec. IIIIB II 

and refers to relativistic corrections, whereas SU^^'' (Sec. lIII B 2t is a 
field-configuration dependent shift. The shift Su^™^^ (Secs. HllB 41 
and llll B 5l refers to radiative corrections to the transition dipole ma- 
trix element, and SU^:!^^ is a dynamic correction to the Rabi frequency 

fSecs. lIIIB3l and lIIIB5t . Finally, is a shift due to corrections 

to the secular approximation (Sec. lIII B6l . 

which has to be much smaller than unity [Xji ^ 1] ■ For hi/2 — 
1000, one obtains 

I(/ii/, = 1000) = 356kHz, (4.14) 
Jnihi/2 = 1000) = 3.6 • 10"*^ , (4.15) 

which means that the probability of one-photon ionization 
does not restrict the above measurement scheme. 

For hi/2 — 1000 and A ~ SOFy^, the theoretical prediction 
for the shift of the Mollow sidebands relative to the central 
Mollow peak by the generalized corrected Rabi frequency is 
as follows: 

±n''^^^^ = ± 100.572258(60) • 10^ Hz . (4.16) 

This formula has been evaluated using the summation formula 
Eq. ( 14. 8> and includes all corrections, in particular the C-term 
evaluated in Sec. IIII B31 For comparison, we also give here a 
theoretical prediction that would be obtained by ignoring the 
C-term, 

±n';^!^c = ± 100.572377(27) • 10^ Hz . (4.17) 



This result is obtained by explicitly setting Cj in Eq. J4.10bl l 
equal to zero, but still using the full summation according to 
Eq. ( 14.91 1 for all other corrections. A comparison of Eq. ( 14.161 1 
to Eq. ( I4.17t shows that the theoretical uncertainties of the 
two results do not overlap. Therefore, the current status of the 
theory would allow to discern the presence or absence of the 
C-term corrections by means of an accurate experiment. 

The principal uncertainty of the generalized Rabi frequency 
in Eq. ( I4.16t is due to the uncertainty in the estimate of the 
nonlogarithmic contribution to Cj and Aj . As a reference, the 
bare shift without corrections is given by 

± Vr22 + A2 ^ ± 99.833975 • 10^ Hz . (4.18) 

This result has to be compared to the radiatively corrected re- 
sult ( I4.16> . 

In Table|I| the individual shifts due to the considered correc- 
tions are listed together with their respective uncertainties for 
both transitions lS'1/2 2P1/2 and lS'1/2 ^ 2P3/2. All shifts 
are evaluated in their "summed" form [see e.g. Eq. ( 13. OH . 

(?) 

For the bare Lamb shift corrections i^ara' uncertainty 
in the shift is mainly due to the uncertainty in the numeri- 
cal value for the Lamb shift of the hydrogen IS" state, see 
Eq. ( I3.10> . The Bloch-Siegert shift acquires a numerical un- 
certainty due to neglected terms of relative order O{fl/oji,) 
from the higher-order Bloch-Siegert-type shifts (higher-order 
perturbation theory in the counter-rotating terms). There is a 
further source of uncertainty for the Bloch-Siegert shifts due 
to terms of relative order ©(A/wl, ^^/^l) from the expan- 
sion leading from Eq. ( 13.191 1 to Eq. ( 13.201 1. The shift due 
to off -resonant excitation T) has an uncertainty because con- 
tributions to the energies Eq. ( I3.26> of relative order A/ljr 
and of order fi/ci^R have been neglected in evaluating the 
t-parameters in Eq. (|^^}. There are also uncertainties of 
the Rabi-frequency shifts due to relativistic corrections to the 
transition dipole matrix elements (Sj), which are due to ne- 
glected higher-order corrections of relative order (Za)'^ with 
respect to the leading corrections. The field-dependent cor- 
rections have an uncertainty due to higher-order effects of rel- 
ative order (Za)^. The main uncertainty of the C-term and 
^j-term corrections are due to the uncertainty which we as- 
sign to the nonlogarithmic contribution of relative parametric 
order a{Za)^ [see Eqs. ( I3.86> ,( l3^92l l. The main uncertainty 
of the shift due to corrections to the secular approximation (S) 
are due to higher-order terms of the expansion leading from 
Eq. ( Il96t to Eq. ( Il97b or relative order 0{A'^/n'^), and of 
fourth-order corrections to the secular approximation of rela- 
tive order 0{T^/Q'^). For the entiies of TableU the shifts have 
been evaluated individually according to Eqs. ( 13.61 1, ( 13.131 1. 
(|I2T), ( Il49l . ( l336l . (im, ( Il72l i. ( ifssl . (E^, and ( 11981 . 

It is perhaps worthwhile to note that one cannot simply add 
the corrections to the quasi-energy of dressed states in the 
same sense as corrections to the energy of bare atomic states. 
For the evaluation of a theoretical Lamb-shift prediction of 
a bare atomic state, the usual procedure is to list the various 
corrections and to simply add these in order to give a theoret- 
ical prediction for the total energy shift (see e.g. Tabs. I and II 
of 1 34]). For the Lamb shift of laser-dressed states, the natural 
interpretation of the corrections implies modifications of the 
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decisive physical parameters that determine the quasi-energy 
of the dressed levels: these are the Rabi frequency and the 
detuning A. This interpretation implies, however, summa- 
tions of the expressions which agree well with the first-order 
results, so that at least for numerically small corrections, it 
is still permissible to simply add the correction terms in an 
approximative sense. The correspondence of first-order ex- 
pressions and summed results holds approximately unless the 
correction is large compared to fl and A; in this case the sum- 
mation yields a different result as compared to the first-order 
expression. This is the case for the numerically dominant ef- 
fect referred to a s SU^^^ ™^^ in Tab.|l] As already discussed in 
Sees. 1111 A 2l and llll A 31 this summation is somewhat nontriv- 
ial in particular for the Lamb shift corrections. Because funda- 
mental symmetry properties prevent the radiative corrections 
from coupling S and P states (even in higher order), prefer- 
ence is given to the summed results. In principle, more ex- 
plicit higher-order calculations outlined in Sec. 1111 A"3] would 
be desirable to verify the summations for all radiative effects 
discussed here. 
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/2 



2P< 



3/2 



In this section, we repeat the above numerical analysis with 
\lSi/2,m=zL'^/2) as the ground state and I2P3/2, m=±Y2) as 
the excited state. The decay width of the 2P3/2 state is given 
by T3/2 = 99.70942(1) • 10^ Hz In order to account for 
the dependence on the laser field intensity, we introduce the 
parameter ^13/2 = |il|/r3/2. Then for ^13/2 = 1000 and A = 
50r3/2, one has tjR S> S> r3/2, A such that the corrections 



to the detuning and the Rabi frequency in Eqs. (|^TOaJ and 
(I4.10b> become 



(V2) 

rad 

(V2) 

rad 



-8.149896(33)- lO^Hz, (4.19a) 
-20.76(56) • 10"^ (4.19b) 



Here, the parameter /13/2 may be expressed in terms of the 
electric field strength £l as 



h,/2 = 490.425 • IQ-^ \£sw [V/m]| 



(4.20) 



For this transition, the ionization to the continuum is again 
given by 

I(/i3/2 = 1000) = 356 kHz, (4.21) 
Jn{h3/2 = 1000) = 3.6 • 10"^ . (4.22) 

Thus the scheme is not restricted by ionization on this transi- 
tion either. 

With the above parameters for the Rabi frequency and the 
detuning, the positions of the Mollow sidebands relative to 
the center component with full corrections, without C-term 
corrections and without any corrections are given by 



^'/'^ = ±100.568846(60) • 10^ Hz , (4.23) 
J^'/'J = ±100.568966(27) • 10^ Hz , (4.24) 
±V'n2 + A2 = ±99.833975 • 10^ Hz . (4.25) 



±17, 



Thus also in this case the theoretical uncertainties of the gener- 
alized corrected Rabi frequency in Eq. j4.23> and of the corre- 
sponding value in Eq. ( I4.24t obtained by ignoring the C-term 
correction shift do not overlap. 

The individual shifts with their uncertainties are listed in 
TableHlas discussed in Sec. lIV Al 



V. DISCUSSION AND SUMMARY 

In this article, we have discussed radiative corrections to the 
usual quantum optical expression for the Mollow spectrum, 
i.e. the resonance fluorescence spectrum of an atom with two 
relevant energy levels driven by a strong coherent laser field. 
To lowest order, the Mollow spectrum consists of one main 
peak which is centered at the frequency of the driving laser 
field and of two sidebands, which are shifted from the center 
by the generalized Rabi frequency fipi = vTP + A^. For the 
analysis, we have used concepts introduced originally in two 
different fields: the dressed-state formalism of quantum optics 
and the renormalized radiative corrections which are treated in 
the formalism of quantum electrodynamics. Throughout the 
analysis, we have adopted the dressed-state basis as the natu- 
ral starting point for our analysis of corrections to the quasi- 
energies of the combined atom-laser system. 

From the point of view of spectroscopy, the Mollow spec- 
trum is attractive for several reasons. First of all, the radia- 
tive corrections manifest themselves in the shift of the Mol- 
low sidebands with respect to the central Mollow peak. Thus 
the tiny radiative corrections are measured relative to the gen- 
eralized Rabi frequency, which for typical parameters of the 
driving laser field is several orders of magnitude smaller than 
optical frequencies. Also, the Mollow spectrum is centered 
around the frequency of the driving laser field. Thus it is a 
kind of a differential spectrum because the laser field detun- 
ing is automatically subtracted. 

As for quantum electrodynamics, radiative corrections to 
the Mollow spectrum are a quantum-field theoretic problem 
in the presence of two classical background fields, the laser 
field and the binding coulomb field. Our analysis also dif- 
fers from typical QED calculations relying on the 5'-matrix 
formalism, as the process under study is dynamical. In order 
to account for the quantum fluctuations of the dipole moment 
leading to the incoherent part of the Mollow spectrum, a static 
description is not sufficient. 

The shifts of the Mollow sidebands may be interpreted as 
arising from corrections to either the detuning or the Rabi fre- 
quency. In particular, the detuning is modified by the bare 
Lamb shift, Bloch-Siegert shifts, and virtual off -resonant ex- 
citations. The Rabi frequency is corrected by relativistic and 
radiative corrections to the transition dipole moment, by field- 
configuration dependent corrections, by a dynamic correction, 
and by corrections to the secular approximation. Of particu- 
lar interest is the dynamical correction to the Rabi frequency. 
This correction arises from an evaluation of the second-order 
radiative self-energy corrections of the combined system of 
atom and laser field in terms of the dressed states of this sys- 
tem. To lowest order of the limit il, A ^ ujeg, this yields a 
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corrections which can be identified with the usual Lamb shift 
of the atomic bare states. The dynamic correction is then ob- 
tained by keeping terms linear in fl, A in the above analysis 
and cannot be explained in terms of the bare state Lamb shift 
alone. 

The corrections to the detuning may be incorporated into 
the analysis by the replacement {j = 1/2, 2/2) 

A^A-A«, (5.1) 

where a|.^|j is defined in Eq. (I4.10a> . Correspondingly, the 
corrections to the Rabi frequency are given by 

+ (5.2) 

The dimensionless quantity Cl'fJ^ is defined in Eq. ( I4.10b> . 
Then, the generalized Rabi frequency supplemented by the 
discussed relativistic and radiative corrections is [Eq. (14. 9H 

In a numerical analysis, we provide a theoretical analysis 
which is required in order to accurately resolve the dynami- 
cal shift. For this, we suppose the driving laser fields to be in 
a standing-wave configuration. As a promising candidate for 
the experiment, we identify the hydrogen 15i/2 2Pi/2 and 
15i/2 ^ 2P3/2 transition. The results are discussed for a driv- 
ing laser field parameter set which is expected to be within 
reach of improvements of the currently available Lyman-a 
laser sources in the next few years. 

For the 151/2 ^ 2P1/2 transition and for fl — 1000 • Ti/^, 
A = 50 • ri/2 , the Rabi frequency is shifted with respect to 
51r = vTP + A^ by relativistic and radiative corrections as 
follows, 

± (n^c'^^^ - r^R.) = ±738.282(60) • 10^ Hz . (5.4) 
The corresponding result for the lS'1/2 ^ 2P3/2 transition with 

fi = 1000 • r3/2, A = 50 • r3/2 is 

± (rj[f/'' - n^j = ±734.871(60) • 10^ Hz . (5.5) 

We note however that we are only concerned with theoretical 
issues. Thus uncertainties due to possible experimental issues 
such as a misalignment of the apparatus or due to additional 
trapping potentials have not been considered. 

In summary, we have presented a detailed analysis of the 
leading nonrelativistic and relativistic corrections to the Mol- 
low spectrum. The analysis includes the relativistic and non- 
relativistic corrections up to relative orders of {Za)"^ and 
a{Za)'^, respectively, and also includes Bloch-Siegert shifts, 
stimulated radiative corrections involving off -resonant virtual 
states, field-configuration dependent corrections and correc- 
tions to the secular approximation. Based on these results, we 
provide a numerical analysis of the corrections of the MoUow 
spectrum of the hydrogen IS — 2P transition. By a compar- 
ison with experimental data, one may verify the presence of 



dynamical leading-logarithmic correction to the dressed-state 
radiative shift, which cannot be explained in terms of the bare 
Lamb shift (see Sec. llIlB"3l . This allows to address questions 
related to the physical reality of the dressed states. On the 
other hand, the comparison with experimental results could 
also be used to interpret the nature of the evaluated radia- 
tive corrections in the sense of the summation formulas which 
lead to the interpretation of the shifts as arising from rela- 
tivistic and radiative corrections to the detuning and the Rabi 
frequency. 
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APPENDIX A: DIPOLE MOMENTS AND SPIN 

The spontaneous emission decay rate of the population of 
an excited state |i) to a final state |/) is given by 

T^Lo^\{f\x\i)\\ (Al) 

where all elements of the position vector have to be considered 
in the coupling with the vacuum field. For the 2P IS de- 
cay in atomic hydrogen, and in the nonrelativistic Schrodinger 
theory without spin, the squared modulus of the dipole mo- 
ment vector is given by 

E 1(1^1 -N2P)P^^(^, (A2) 

independent of the magnetic quantum number 
nii E {—1,0,1} of the "initial" P state (by the term 
"final state" we will denote in this section the particular state 
that enters as a "bra-" in the Dirac notation, i.e. the IS* state 
in the above case). In the Schrodinger-Pauli theory with spin, 
one has different decay channels depending on the spin state 
of the initial and the final state. For the channel where the 
initial and the final state have the same magnetic quantum 
number nii — nif — one obtains 

E 1(1%., m.| |2Pv.,m.)p = ^ ^ -J^ . (A3) 

In contrast, the channel with opposite magnetic quantum num- 
ber mi = — m/ = ±1/2 yields a contribution of 

g |(15v2,-™.| |2Pv.-.)p = - ^ . 

(A4) 
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This calculation predicts that if one were to measure the elec- 
tron spin polarization in the final state, then the \2Pi/^ state 
would be twice as likely to decay into a llS'1/2) state with 
opposite total electron angular momentum than into a llS'1/2) 
state with the same total electron angular momentum as the 
initial \2Pij^) state. Adding the two decay channels to the fi- 
nal state, we obtain 

3 

1 



3I" {Zafm?- ' 



(A5) 



i.e. the same result as in the spinless case, as it should be. For 
the decay of the 1 2P3/2 , = ± 1/2) state, one obtains the same 
total decay rate, but here the decay with to/ = mj is twice as 
likely as the decay with to/ = —to^. 

These results for the dipole moments have to be reconsid- 
ered for excitation of an atom in the llSi/a, m^i) ground state 
with a laser field which is linearly polarized in one direction, 
say the z-direction. Then other than for the interaction with 
the vacuum field, only the z-component of the dipole moment 
vector has to be considered. These matrix elements can be 
used to calculate the Rabi flopping frequency corresponding 
to the driving laser field. For the spinless case, one obtains 



(A6) 



We now include the spin and choose a definite initial state 
\lSi/2,rai — +1/2). We obtain 



2 1 



3 Zam 



(A7) 



where we have omitted the "m^ ("to / =") from the initial 
(final) state vector Furthermore, one has 

(151/2, V2U |2Pv2,-V2) =0, (A8) 
(lSi/2,V2|a:|2Pi/2,V2) =0, (A9) 
(lSi/2,V2|y|2Pi/2,V2) = 0, (MO) 

so that the z-polarized field only couples the 2P1/2 state with 
TOi = TO/ to the ground state, and this excited state is only 
one coupled to the ground state by z-polarized light. For the 
I2P3/2, TO/) upper state, the corresponding results are 
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(lSl/2,V2|.|2P3/2,V2) = |,V3Z^' ^^^^^ 

(15i/2,V2|^|2P3/2,-V2) -0, (A12) 

(15i/2,V2|2;|2P3/2,V2) =0, (A13) 

(lSi/2,V2|2; |2P3/2,V2) =0, (A14) 

(lSi/2,V2|2|2P3/2,±3/2) =0. (A15) 

Thus also in this case only the upper state with same mag- 
netic quantum number is coupled to the ground state, but with 



a matrix element which differs by a factor of — \/2 in mag- 
nitude from the corresponding result for the |2Pi/2, 1/2) upper 
state. From these spin-resolved results, the corresponding ma- 
trix element without spin Eq. ( IA6> may be obtained by sum- 
ming over the final states and averaging over the initial states. 



APPENDIX B: EVALUATION OF THE MATRIX ELEMENTS 

In this section we demonstrate the evaluation of the matrix 
element in Eq. (|^2}: 



G {O . 



(Bl) 



where |e) is the 2P, to = state. We start by calculating the 
"unreduced" matrix element 



Me(C) = (e|zG(C)z|e) 



(B2) 



where the full sum over intermediate states is employed in 
G(C), and the wavefunction of the excited state state $2_p,m=o 
is given by a product of a radial and an angular contribution: 

$2P,m=0(r) =i?2pWyio(0,</'). 

Here, (r, 9, (p) are spherical coordinates. In these coordinates 
a representation of the Green function in position space is 
given by [45^, Eq. 2.2] 



where g„ (ri ,r2,i^) is the radial component of the 
Schrodinger-Coulomb propagator 



2m 



21+1 



(rir2)' exp 



ri + 7-2 



E 

fe=0 



fc! i2i+l / 2r^\ ^2/+l ( 2r^' 



{2i + i + ky.{i + i + k-i^) 



(B3) 



containing associated Laguerre polynomials L'^{r). The 
quantity ly = (03%/— 2toE)^i is an energy parameter which 
is related to the parameter t used in Sec. 1111 A5l bv v = nt 
where n is the principal quantum number of the initial bound 
state. The Bohr radius is defined in ( I3.44t . and we evaluate all 
matrix elements here for the case Z = 1 (atomic hydrogen). 
Thus for the 2P state discussed here we have v — 2t. The 
index / is summed over all possible angular momentum num- 
bers of the virtual intermediate states. Starting from a P state 
with I — 1, both S (I — Q) and D (I — 2) states are possible as 
intermediate states. The integration may be further separated 
in angular and radial parts: 



1=2 , 



where the angular integrations yield 



M 



1=0 



M, 



=2 



4 
15 



(B4) 



(B5) 
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The radial parts may written explicitly as 

/•OO 

/ dridr2rlrlR;p{ri)R2p{r2) 
Jo 



Ml = 



X gi{ri,r2,v) , 
for Z = 0, 2. Simplifying further, one obtains 

fc! 



(B6) 



12ta? 



fe=0 



(l + fc)!(l + A:-2i) 



dr r L 



\aBtJ 



(B8) 



The integral in /o can be evaluated using (see 
Sec. 6.10] and |52, Sec. 10.12]) 



^ ^ ^ 7 TM/ A A-l-Tr(7+l) 



n!r(/x+ 1) 

X r(^ + n + 1) 2Fi(-n,7+ 1,^+ 1, A^i) 



to give 



(l+t)s 



fc! 



1 + i 



(B9) 



(BIO) 



Using an explicit expression for the hypergeometric func- 
tion IdU Sec. 2.1.1], we obtain 



2Fi{a,b,c;z) ^ ^ 



3=0 



(BID 



where the Pochhammer symbols {a)j are given by 

r(a + j) 



T{a) 



(B12) 



Contiguous relations for the hypergeometric function jsTl 
Sec. 2.8] then lead to 



M 



;ad 



14J-11 



2'H 



Xo{t) = 45 - 90i - 84i2 + 258^"^ + 18t^ - 294t^ 
+148i^ - 2f + 257t^ , 



where the hypergeometric function <l>(n, is defined in 
Sec. 1111 A 51 A similar calculation for the intermediate D 
{I = 2) states yields 



M 



■rad 
1=2 



16 t^Mt) 

3(t-l)7(t + l)5 

A'2(i) = -45 + 90t + 165 f - 420<3 - 174i* + 768t^ 
-3At^ - 700f - - 121 At^ + 4733i^° . 



Inserting in Eq. ( IB4t finally yields the expression in Eq. 
J3.41bl l. The reduced matrix element can then be obtained 
from this by subtracting the contributions of the two interme- 
diate states |e), \g). The excited state contribution vanishes 
due to parity, and the ground state contribution is given by 



l(gk|e)|^ 

c 



E- 



(B13) 



IS 



This term, which cancels the divergence as C — > Eis in ( I3.46> . 
may be verified by inserting the resonant 1 5 state as the inter- 
mediate state into the matrix element. Alternatively, the can- 
cellation may be seen as follows: On setting the intermediate 
state energy to Eis + e, the series expansion of the unreduced 
matrix element Me{Eis) receives a contribution proportional 
to 1/e which diverges for e ~> 0. This diverging part is can- 
celed by the intermediate state contribution Eq. jB13l l in the 
reduced matrix element A/g {() to give a finite result. 

If one compares this derivation with a similar calculation 
for the standard matrix element 



M = ^(e|x*G"(C)a;' 



(B14) 



rang -;— rrad rang -;— rrad 

= Ml Jo ■ Mi^o + M1J2 ■ Mi^2 , (B15) 



where all polarization directions are considered, one finds that 
the respective radial parts for I = and I = 2 are identical to 
the ones in Eq. (IB4> : 



. . rad ^r-rrad 

Mi=o = Mi=o, 



-T— jrad -;p-7rad 



(B16) 



For the angular parts however one finds 



M 



1=0 



M' 



1=0 ' 



M 



1=2 



5 ^r-ran, 

- M 
2 



1=2 ■ 



(B17) 



The reason for this is that the fixed polarization of the driv- 
ing laser field only allows to excite one of the magnetic sub- 
levels of the intermediate S and D states. The sum over all 
polarizations in the standard matrix element still only gives 
one magnetic sublevel for the intermediate S states, but three 
possible virtual D states. Due to this asymmetry the desired 
matrix element Eq. ( I3.41b> cannot be calculated directly from 
the standard matrix element in Eq. ( IB14t . 
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